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ABSTRACT 

The combination of STIS and optical spectroscopy with STIS imaging of the field surrounding 
J2233— 606 afford a unique opportunity to study the physical nature of the Quasar Absorption Line 
systems. We present an analysis of the ionization state, chemical abundances, and kinematic charac- 
teristics of two metal-line systems at z = 1.92 and z — 1.94 toward J2233— 606. We focus on these 
two systems because (i) the observations provide full coverage of the Lyman series, hence an accurate 
determination of their HI column densities and (ii) they exhibit many metal-line transitions which allow 
for a measurement of their ionization state and chemical abundances. 

Line-profile fits of the Lyman series for the two systems indicate log AT (HI) = 17.15 ±0.02 cm -2 evenly 
distributed between two components for the z = 1.92 system and log A(HI) = 16.33 ± 0.04 cm -2 for the 
z = 1.94 system. By comparing observed ionic ratios of C and Si against calculations performed with the 
CLOUDY software package, we find the ionization state is high and well constrained in both systems. 
Applying ionization corrections to the measured ionic column densities, we determine that the systems 
exhibit significantly different metallicities: w 1/50 solar and < 1/200 solar for the two components of 
z = 1.92 and rs 40% solar at z = 1.94. The properties of the z = 1.92 system are consistent with a 
low metallicity galaxy (e.g. a dwarf galaxy) as well as absorption by large scale structure of the IGM. 
On the other hand, the high metallicity of the z — 1.94 system suggests significant star formation and 
is therefore most likely associated with a galactic system. Furthermore, we predict it is the most likely 
to be observed with the STIS imaging and follow-up observations. If a galaxy is identified, our results 
provide direct measurements on the properties of the ISM for az«2 galaxy. 

Subject headings: quasars: individual (J2233— 606) — absorption lines — galaxies: abundances 



1. INTRODUCTION 



Ide ntifying the physical nature of the Quaoai 1 Aboorp 

tion Lino — (QAL) — oyotoma will have oignificant impact 

on our understanding of galaxy formation, chemical evo- 
lution, and large scale structure (LSS). At high red- 
shift, where QAL systems offer one of the most efficient 
means of probing the early universe, any insight is no- 
table.i To really utilize QAL samples as tools to study 



cosmo logy (LSS, dN/dz) and early galaxy formation, wo 

need to relate the absorbers to astrophysical environ- 
ments and objects. With the advancement of modern tele- 
scopes and instrumentation, in particular high resolution 
echelle spectroscopy, observers are now accurately exam- 
ining the physical characteristics of th ese systems. For th e 
damped Lya systems, the metallicity (Pcttini ct al. 1997), 



1 Based in part on observations with the NASA/ESA Hubble Space 
Telescope obtained at the Space Telescope Science Institute, which 
is oper ated by AURA) Inc. undor NASA Contract NAS 5 26555j 



and ot scrvations collected at European Southern Observatory, La 
Silla, Chile (ESO Nr. 60.B-0381); and observations collected at the 
Anglo- Australian Observatory. 



abun dance patterns (Lu et al. 1996, Prochaska fc Wolfe 



1999 ) and kinematic characteris tics ( Prochaska fc Wolfe 
1997b , Prochaska fc Wolfe 1998 ) have been well studied 
and have revealed vital clues to their physical nature. 

Recent observations of the Hubble Deep Field South 
afford a unique opportunity to investigate the physical 
nature of the QAL systems. With the combination of 



STIS spectroscopy and S TIS imaging ( Gardner ct al, 
1999j |Williams ct al. 1999| of the fiel d toward J2233-606 



(z em — 2.238, B — 17.5; Boyle 1997 ) researchers may di 
rectly identify the absorbers responsible for the QAL sys- 
tems. In turn, these observations will further develop the 
relationship between QAL systems and galaxies or large 
scale structure at moderate redshift. In this paper, we 
investigate the physical properties of the two metal line 
systems at z w 1.9 toward J2233— 606. We choose these 
two systems because the STIS spectroscop y coupled with 



optical echelle spect roscopy of J2233— 606 (Outram et al 



1998 , Bavaglio 1998 ) allows one to well constrain the ion- 
ization state of these systems. By applying ionization cor- 
rections to the measured ionic column densities, we present 
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Table 1 
Spectroscopic Data 



Instrument /Telescope 


Coverage 


Resolution 


Source 


E230M-STIS/HST 


2300- 


-3100 A 


R = 30000 


HST HDF-S STIS Team 


G430M-STIS/HST 


2300- 


-3100 A 


0.56 a 


HST HDF-S STIS Team 


UCLES/AAT 


3530- 


-4390 A 


R = 35000 


Outram et al. (1998) 


EMMI/NNT 


4386- 


-8270 A 


14 km/s 


Savaglio (1998) 



chemical abundances for Si, C, N, and Fe. Finally, we 
comment on the kinematic characteristics of the absorp- 
tion line systems and speculate on the systems' association 
with protogalaxies and large scale structure at high red- 
shift. 

In § 2, we discuss the multiple data sets used in this 
analysis. Section 3 details our fits to the Lyman series of 
the two QAL systems. The ionic column densities, ioniza- 
tion state, and chemical abundances are examined in § 4. 
Finally, we speculate on the impact of these measurements 
in § 5. 

2. DATA 

In the following analysis we utilize several data sets ob- 
tained with a number of instruments and telescopes (Ta- 
ble |]). These include ultra-violet spectra taken with the 
STIS E230M echelle and G430M grating on board the 
Hubble Space Telescope, and optical echelle spectroscopy 
acquired with the AAT and NTT and kindly provided 
by Outram et al. (1998) and Savaglio (1998). In all, 
this collection of data provides continuous coverage from 
2275-8300A with resolution R « 30000 (except the region 
from 3100 — 3500A corresponding to the G430M grating) 
and moderate signal to noise (S/N). With the exception of 
the NTT spectrum where we have adopted the continuum 
provided by Savaglio, we have continuum fit the spectra 
with low order Legendre polynomials using a package sim- 
ilar to the Iraf routine continuum. 



3. HI MEASUREMENTS 

A requisite facet in analyzing the ionization state and 
chemical abundances of QAL systems is an accurate deter- 
mination of the HI column density. For the QAL systems 
at z = 1.92 and z = 1.94 towards J2233-606, the echelle 
observations performed with STIS on board the HST are 
essential as they allow for a fit to the higher order lines 
in the Lyman series. We have measured the HI column 
densities by fitting Voigt profiles to the Lyman series with 
the VPFIT package kindly provided by R. Carswell and 
J. Webb. Figures [y and || and Table || present the results 
of this analysis. Note that in performing the fit to the 
Lyman series, we have included fits to other Lya forest 
clouds when necessary. 

In the z=1.92 system, most of the H I is split between 
two main components at z =1.9256290 and 1.9279202. 



The total column density is logiV(HI) = 17.15±0.02 cm -2 
with a nearly even split between the main components. 
As suggested by Outram et al. (1998), the STIS observa- 




-200 200 400 

Relative Velocity (km s _1 ) 

Fig. 1. — A plot of the Lyman series for the Lyman limit 
system at z = 1.92 toward J2233— 606. Overplotted is the 
VPFIT solution which includes several coincident Lya for- 
est clouds. The dotted regions indicate data not included 
in our fit to the Lyman series. 
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tions clearly demonstrate that the Lyman break evident 
at « 2700 A is due to the Lyman limit (LL) system at 
z = 1.92. For the z=1.94 system, the majority of H I falls 
at z=l. 9425817 and the column density is well constrained 
by the high-order unsaturated Lyman lines (Ly-10, 11, 12). 
We find logA(HI) = 16.33 ± 0.04cm- 2 . Therefore this 
system has a minimal contribution to the observed Lyman 
break. 



Table 2 

hi fits for the z=1.92 and z=1.94 qal systems 




-200 -100 100 200 

Relative Velocity (km s _1 ) 



Fig. 2. — A plot of the Lyman series for the QAL system 
at z = 1.94 toward J2233-606. Overplotted is the VP- 
FIT solution which includes several coincident Lya forest 
clouds. The dotted regions indicate data not included in 
our fit to the Lyman series. 



Only weak limits can be placed on D/H in these ab- 
sorption systems. In the z=1.94 system, extra H I falls 
to the blue of the main H I component (at —60 km/s) 
making future detection of the D I line impossible. In the 
z=1.92 system, the Lya line is too wide for a detection or 
an upper limit. The blue wing of Ly/3 might be used for 
an upper limit, but only a weak limit is derived from the 
STIS E230M spectrum: D/H < 10" 3 . 



z 


log N 




b 


°~b 




( cm~^ 

V I 


( cm - ^ 


(km s — * ) 


(km s — * ) 


1.2574902 


12.93 


0.17 


8.1 


4.0 


1.2847255 


13.81 


0.05 


48.5 


6.4 


1.3393718 


13.28 


0.09 


42.6 


10.2 


1.3529782 


13.42 


0.07 


21.1 


3.8 


1.4667628 


13.77 


0.06 


24.9 


2.5 


1.9407216 


13.80 


0.09 


49.2 


6.1 


1.9419496 


15.46 


0.06 


20.3 


1.6 


1.9425205 


16.33 


0.04 


22.4 


1.0 


1.9433281 


13.81 


0.08 


37.4 


8.7 


1.9256160 


16.79 


0.03 


33.3 


1.2 


1.9269993 


15.69 


0.05 


44.7 


4.9 


1.9277830 


16.76 


0.06 


21.7 


1.8 


1.9281701 


16.09 


0.19 


23.0 


4.8 


1.9288300 


15.37 


0.05 


25.1 


1.1 



4. IONIZATION STATE AND ABUNDANCES 

Since the discovery of QAL metal-line systems, re- 
searchers have investigated the ionization state of these 
systems by comparing observations of multiple ions for 
a given element against theoretical predictions from pho- 
toionization models (e.g. Bteidcl 1990, Hamman fc Ferlanc| 
1992). While the initial analyses placed meaningful con- 



straints on the ionization state, they were limited by less 
accurate column density measurements derived from lower 
resolution data. Recently, Prochaska (1999) has demon- 
strated that one can precisely determine the ionization 
state of a LL system given the very accurate column den- 
sity measurements afforded by modern telescopes. The 
combined HST and optical echelle spectra for J2233— 606 
provide the opportunity for such an analysis. In what fol- 
lows, we investigate the ionization state of the two metal- 
line systems at z ~ 1.9 toward J2233— 606 under the as- 
sumption that they are photoionized by the extragalactic 
ultra-violet background (EUVB) radiation. We adopt two 
input spectra: (i) the Haardt-Madau (1996; HM) EUVB 
spectrum which consists entirely of radiation from back- 
ground quasars attenuated by the IGM and (ii) a model 
compr ised of a synthesis of backgrou nd galactic and quasar 
light QValls-Gabaud fc Vcrnct 19991 ; VV) also attenuated 
by the IGM. The latter model predicts a softer spectrum 
owing to the absorption of radiation with hv > 4 Ryd by 
the ISM of the emitting galaxies. 

Central to our analysis is the CLOUDY software pack- 
age (V 90.04) developed and kindly provided by G. Fer- 
land (1995). We assume a plane-parallel geometry and 
approximate an isotropic background by placing a point 
source at a very large distance. There is a homologous 
relation between the intensity of the EUVB radiation and 
the volume density of the QAL system which is typically 
parameterized via the ionization parameter, U . Here, we 
define the ionization parameter by relating the intensity 
of the EUVB at 1 Ryd (J 912 ) with n H - 
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U = *™ = = (2 x 10- 5 ) J «"/ 10 " 2 '- 5 (i) 

cnn AnhcriH n#/cm 6 

The CLOUDY software package inputs the photoioniza- 
tion spectrum, U, and the observed A(HI) and then cal- 
culates the mean ionization fraction for all relevant ele- 
ments. This enables one to precisely determine the ion- 
ization state of a QAL system and then apply ionization 
corrections to the observed ionic column densities to de- 
rive chemical abundances. We turn now to investigate the 
QAL systems individually. 



0.5 




0.5 



Felll 1122 
1 I 1 1 1 1 l~ 




SUV 1393 



0.5 - 



0.5 



0.5 



m 



Sill 1260 



CII 1334 





Relative Velocity (km s ) 

Fig. 3. — Metal-line profiles for the LL system at z — 1.92. 
Overplotted are our fits to the absorption profiles derived 
with the VPFIT software package. For completeness, we 
include all metal-line transitions where one might expect 
to observe an absorption profile. We indicate regions of the 
profiles contaminated by other absorption lines by dotting 
the data. The error array is the light gray line. 



4.1. z=1.92 

This QAL system has a total HI column density iV(HI) = 
10 1715 cm -2 and is responsible for the Lyman limit at 
« 2700 A. It exhibits a number of metal-line transitions 
commonly observed in LL systems (presented in Figure ||) . 
To determine the ionic column densities from the absorp- 
tion line profiles, we have utilized the VPFIT software 
package in a manner similar to that for the iV(HI) deter- 
minations. When fitting the profiles we introduced the 
minimum number of components necessary to yield a sta- 
tistically 'good' fit. The profiles trace one another very 
closely in velocity space such that we successfully tied the 
redshifts and Doppler parameter of each component for 
all of the ions. This point has significant bearing on our 
analysis because we believe the correspondence indicates 
all of the ions are spatially coincident and can be treated 
together in a single photoionization calculation. Further- 
more, it implies there is little variation in the ionization 
state along the sightline. The VPFIT solutions and 1<7 
errors are presented in Table |[ With only a few excep- 
tions, the solutions very closely match the results from 
Outram et al. (1998). We have also checked the ionic 
column density measur ements by applying the Ap parent 
Optical Depth Method (Savage and Sembach 1991) which 
shows excellent agreement. For those transitions which ex- 
hibit no significant absorption, we report 3a upper limits 
to their ionic column densities. 



Table 3 
Ionic Column Densities for the 



= 1.92 System 



Comp 


z 




b 




Ion 


log N 


°"logJV 






icr 6 


km/s 


km/s 




cm 


cm -2 


1 


1.925970 


3 


6.13 


0.41 


Si+ 3 


12.78 


0.02 








6.13 




Si++ 


13.84 


0.78 








9.37 




c+ 


13.12 


0.03 








9.37 




C +3 


13.12 


0.14 








8.68 




N++ 


13.79 


0.12 








8.68 




N+ 4 


13.06 


0.15 


2 


1.925695 


14 


13.89 


3.00 


Si+ 3 


12.25 


0.08 








13.89 




Si++ 


12.32 


0.38 








21.24 




c+ 


12.78 


0.10 








19.67 




N++ 


12.93 


0.62 












Fe+ 


< 13.3 














Fe++ 


< 13.5 





The fact that the HI gas is distributed between two main 
components complicates the analysis. Examining Figure [}], 
one notes that the metal-line profiles lie near the blue HI 
component (z — 1.9256), albeit not strictly at the center 
of the profile (compare with the absorption profiles for the 
z = 1.94 system). Furthermore, we observe no significant 
metal-line profile for the HI component at z = 1.9278. 
It is very possible that the LL system is actually com- 
prised of two spatially separated absorbers with a velocity 
separation of w 200 km s -1 along the sightline. Because 
we cannot constrain the geometry of these absorbers, our 
CLOUDY simulations are an oversimplification. In the fol- 
lowing, we treat the absorbers separately and focus on the 



PROCHASKA & BURLES 



5 




Fig. 4. — The results from a series or ^LOUDY calculations (Ferland 1995J' assuming logTV(HI) = 16.8 with (a) a 
Haardt-Madau ionizing spectrum and (b) the Valls-Gabaud/Vernet model. Plotted are the calculated mean ionic column 
densities (assuming intrinsic solar abundances and [Fe/H] = —0.5) versus a modified ionization parameter, U, as defined 
in the text. The vertical dashed lines indicate the observed C + /C +3 ratio (with errors) and the dotted vertical line is a 
lower limit to U determined from the upper limit to the Si + /Si +3 ratio. 



HI component at z = 1.9256 with JV(HI) = 10 1 
which corresponds to the metal-line profiles. Finally, we 
will use the results of our analysis to place limits on the 
metallicity of the other HI component. While it is possible 
that the z = 1.9278 absorber may attenuate the EUVB 
radiation, we find that the ionization corrections hardly 
vary for logA(HI) = 16.8 — 17.2 cm -2 and therefore do 
not consider shielding to be a problem. 

Table 4 

Chemical Abundances for the z=1.92 System 



Ion 


Haardt/Madau 


Valls-Gabaud/Vornot 




[X/H] 


[X/H] 


C+ 


-1.47 ±0.1 


-1.70 ±0.1 


C +3 


-1.47 ±0.2 


-1.70 ±0.2 


N++ 


-1.35 ±0.2 


-1.66 ±0.2 


N +4 


-2.28 ±0.3 


-2.57 ±0.3 


Si+ 


< -1.54 


< -1.68 


Si+ 3 


-1.16 ±0.15 


-1.43 ± 0.15 


Fe+ 


< 1.6 


< 1.5 


Fc++ 


< -0.1 


< -0.2 



Given the ionic column densities and our measurement 
of iV(HI), we can constrain the ionization state of this sys- 
tem by comparing the relative ionic column densities with 
the calculations from the CLOUDY package. To eliminate 



the effects of the underlying abundance pattern (e.g. Type 
II SN enhancements) and dust depletion, we will focus on 
multiple ions from a single element. For this LL system, 
the tightest constraints are imposed by the C + /C +3 ratio. 
Figure || presents the CLOUDY solutions (calculated ionic 
column densities) for a plane-parallel slab with A(HI) = 
10 16 ' 8 cm" 2 photoionized by (a) an HM EUVB spectrum 
and (b) the VV spectrum for a range of ionization param- 
eter assuming [Fc/H]=— 0.5 dex. The results are largely 
independent of the metallicity assumed in the CLOUDY 
calculation; varying [Fe/H] uniformly varies all of the pre- 
dicted column densities. The dashed vertical lines indi- 
cate the U values corresponding to log[JV(CII)/JV(CIV)] = 
0.16 ± 0.2 while the dotted line denotes a lower limit to 
U derived from log[7V(SiH)/7V(SiIV)] < -0.66. For the 
HM spectrum, the C + /C +3 and Si + /Si +3 ratios are con- 
sistent at the 2(7 level and indicate logc7 = —2.69 ± 0.1 
and N(H) = 19.21 ± 0.20 dex. This implies n H = 
1 x 10~ 2 cm~ 3 ( J 31 2/10~ 2 ^ 5 ) which ag r ees with limits t o 
Uh in other LL systems ( Bteidel 1990 , Prochaska 1999 ). 
In terms of the VV EUVB model, the results are com- 
parable with a slightly higher prediction for A(H) and 
U: logU = -2.38 ±0.1, N(R) = 19.51 ± 0.2 dex, and 
nn = 5 x 10 -3 cm -3 (Jgi 2 /10 _21 - 5 ). Ignoring the effects 
of clumping within the LL system, this implies a charac- 
teristic length scale, I ~ N(H)/nu ~ 1 kpc. 

Having placed relatively tight constraints on the ioniza- 
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tion state, we can calculate the chemical abundances for 
this system by making ionization corrections to the ob- 
served ionic column densities. Table | lists the abundances 
relative to solar ([X/H] = log[N(X)/N(H)] - log[X/H] ) 
derived from each ion for the two EUVB models. The er- 
ror bars account for uncertainties in AT(HI), U, and the 
ionic column density. Note the very small error to the 
[C/H] value derived from A(CII). As one observes from 
Figure ||, the relevant U values span a very flat section 
of the CII curve such that [C/H] is nearly independent 
of U. Therefore, the Carbon abundance has a very small 
error to its ionization correction, i.e., the error in [C/H] 
is dominated by the uncertainty in the AT(HI) measure- 
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Fig. 5. — Metal-line profiles for the QAL system at 
z = 1.94. Overplotted are our fits to the absorption pro- 
files derived with the VPFIT software package. For com- 
pleteness, we include all metal-line transitions where one 
might expect to observe an absorption profile. We indicate 
regions of the profiles contaminated by other absorption 
lines by dotting the data. The error array is the light gray 
line. 



ment. This implies CII measurements can provide reason- 
ably robust metallicity measurements for all LL systems 
with AT(HI) w 10 17 cm- 2 . For this system, we find [C/H] 
= -1.47 ± 0.1 dex for the HM spectrum and [C/H] = 
— 1.70 ± 0.1 dex adopting the VV model. The values are 
ty pical of other estimates for Lyman limit systems at high 
z ( |Stcidel 1990| , [Buries fc Tytler 1998j ). Comparing the Sil- 
icon abundance ([Si/HJ = -1.16±0.15 for HM and [Si/H] 
= -1.43 ± 0.15 for the VV model) derived from V(SiIV) 
against [C /H] we observe good agreement which gives fur- 
ther confidence in our modeling as the values are indepen- 
dent of the C and Si ionic ratios. The relative abundance 
of C to Si, [C/Si] = — 0.3±0.1 dex, is consistent with solar 
abundances at the 3tr level, but favors the ratio of —0.3 dex 
observed in Population II stars as expected for low metal- 
licity systems. In all, we confidently report that the HI 
component at z = 1.9256 has a metallicity of « 1/50 solar 
abundance. Without any ionic column density measure- 
ments with which to constrain the ionization state of the 
HI component at z = 1.9278, one can only speculate on 
its metallicity. Assuming that the ionization state matches 
that of the blue HI component, we find [C/H]< —2.4 dex 
from a 3cr upper limit to the CII A1334 column density. 



4.2. 



-1.94 



From an examination of a low resolution spectrum of 
J2233-606, Sealey et al. (1998) proposed this QAL sys- 
tem is responsible for the Lyman limit system observed at 
» 2700 A. Fi gure clearly demonstrates, however, that 
the Lyman series is not saturated at Ly-10 and therefore 
the system is optically thin at the Lyman limit. In Fig- 
ure we present the velocity profiles of the metal-line tran- 
sitions overplotted by the VPFIT solution (Table |^) for 
those transitions where an unblended absorption feature 
is present. As with the Lyman limit system at z = 1.92, 
the low ionization states (e.g. C + , Si + ) very accurately 
trace the high-ion profiles. Therefore, we fit all of the pro- 
files by tying the redshifts and Doppler parameters of the 
individual components. Unlike the system at z = 1.92, the 
metal-lines also closely trace the HI gas and we treat all 
of the gas as a single system. 

Figure || presents the CLOUDY calculations for a QAL 
absorption line system with JV(HI) = 10 16 ' 33 cm -2 and 
[Fe/H] = -0.5 dex adopting (a) a HM EUVB spectrum 
and (b) the VV EUVB model. Since this system is op- 
tically thin at the Lyman limit, the results are largely 
insensitive to the iV(HI) value. Similar to our exam- 
ination of the z — 1.92 system, we overplot the al- 
lowed values for the ionization parameter given the ob- 
served ionic ratios for Si + /Si +3 and C + /C +3 . Because 
the Silll and CIV profiles are heavily saturated, the 
most significant constraints for the ionization state are 
from the Si+/Si+ 3 ratio, V(SiII)/V(SiIV) = -0.70 ± 
0.15 dex. For the HM spectrum, we find logf/ = —2.46 ± 
0.1 which implies logAf(H) = 18.98 ± 0.15 cm- 2 and 
njjRi6x 10 -3 cm~ 3 (J912/IO- 215 ). Meanwhile, compar- 
ing against the CLOUDY solution for the VV EUVB spec- 
trum, the values are: logf/ = -2.24 ± 0.1, logiV(H) = 
19.19±0.2 cm" 2 , and n H « 6 x 10~ 3 cm~ 3 ( J912/IO" 215 ). 
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Fig. 6. — The results from a series ofTJLOUDY calculations (Ferland 1995) assuming logiV(HI) = 16.33 with (a) a 
Haardt-Madau ionizing spectrum and (b) the Valls-Gabaud/Vernet model. Plotted are the calculated mean ionic column 
densities (assuming intrinsic solar abundances and [Fe/H] = —0.5) versus a modified ionization parameter, U, as defined 
in the text. The vertical dotted lines indicate the observed Si + /Si +3 ratio (with errors) and the dotted vertical line is a 
very rough estimate to U from the C + /C +3 ratio. 



Table 5 

Ionic Column Densities for the z=1.94 System 



Comp z a z 


b 




Ion 


log TV cr logJ v 


10" 6 


km/s 


km/s 




cm cm 


1 1.942614 4 


6.17 


0.44 


Si+ 


12.69 0.12 



6.17 
6.17 
9.44 
9 ?1 



Si++ 
Si+ 3 

n+3 



15.03 
13.32 
13.03 
i 5 ^ 



1.942441 



10 



1.941976 



8.74 
11.11 
11.11 
11.11 
17.00 
17.00 
15.76 
11.26 
11.26 
17.22 
17.22 



0.83 



0.53 



N++ 
Si+ 
Si++ 
Si+ 3 

c+ 

C +3 

N++ 

Si+ 

Si+ 3 

c+ 

C+ 3 

N+ 4 
Fe+ 
Fe++ 



13.67 
12.71 
13.06 
13.44 
13.56 
13.80 
13.71 
12.15 
13.03 
13.26 
14.04 

< 13.3 

< 13.1 

< 13.1 



0.81 
0.06 
0.04 
n <^p. 



0.18 
0.10 
0.43 
0.04 
0.05 
0.27 
0.14 
0.20 
0.03 
0.03 
0.14 



tion corrections for Si + and C + are very small (< 0.1 dex) 
so it is very reassuring that the derived [C/H] and [Si/H] 
values are in good agreement: [C/H] = —0.27 ±0.1 (HM), 
[C/H] = -0.50 ±0.1 (VV) and [Si/H = -0.14 ±0.1 (HM), 
[Si/H] = -0.31 ± 0.1 (VV). Altogether, the system has 
a relatively large metallicity, w | — | solar metallicity. 
This value is even larger than what is typi cally observed 
for th e damped Lya systems at z ~ 2 ( Pettini ct al. 
1997| , |Prochaska fc Wolfe 199"9| ). In § 4 we will discuss 



the implications of this result. 

Table 6 

Chemical Abundances for the Z=1.94 System 



Ion 


Haardt /Madau 
[X/H] 


Valls-Gabaud/Vernet 
[X/H] 


c+ 


-0.27 ±0.1 


-0.50 ±0.1 


C +3 


-0.06 ±0.6 


-0.15 ±0.6 


N++ 


-0.93 ±0.2 


-1.17 ±0.2 


N +4 


< -3.00 


< -3.15 


Si+ 


-0.14 ±0.1 


-0.31 ±0.1 


Si+ 3 


-0.14 ±0.1 


-0.31 ± 0.15 



Table B lists the abundances relative to solar for this sys- 
tem for the two EUVB models. The variance in the ioniza- 



4.3. HIGHER IONIZATION GAS 



INVESTIGATING THE METAL LINE SYSTEMS AT z=1.9 TOWARD J2233-606 IN THE HDF-S 



We have also considered constraints on the ionization 
state for each system placed by a lower limit to the 



systems are asso ciated with the outer regions (e.g. halos) 
of w L* gal axies (Bergeron fc Boisse 199lL |Steidel fc Dick 



N ++ / N +4 ratio 2 . The Nitrogen ratios are inconsistent inson 1992 ). This suggests that at z ss 2 one expects the 



with both EUVB models for all values of U in each system. 
We expect this disagreement arises because of a mistreat- 
ment of N+ 4 . N+ 4 is the only ion we examined with an 
ionization potential greater than 4 Ryd and is therefore 
the most sensitive to three factors: (1) the shape of the 
EUVB spectrum at hv > 4 Ryd, (2) the optical depth 
of Hell in the IGM, and (3) the optical depth of Hell 
within the metal-line system. Because the number den- 
sity of quasars peaks near z — 2, it is unlikely that the 
EUVB spectrum is significantly softer than the VV model 
which does account for a contribution from background 
galaxies. Regarding the second point, the optical depth of 
Hell in the IGM is measured to be quite small at z = 2.4 
(jHe.ii — 1-0 ±0.07; Davidson et al. 1996 ) and is expected 
to be significantly lower at z w 2. Therefore, we expect the 
high N ++ /N +4 ratio implies an oversimplification in mod- 
eling Hell within the CLOUDY framework. Whereas the 
CLOUDY package assumes the metals are uniformly dis- 
persed throughout the absorption-line system, it is quite 
possible the metals are primarily located in the central re- 
gions where they could be shielded from high energy pho- 
tons by the surrounding He gas. For both systems, one 
predicts AT(Hell) > 10 17 8 cm~ 2 which imply an optical 
depth THe.il > 1 that is more than enough to account for 
the ss 0.3 dex discrepancy between the observations and 
CLOUDY predictions. Unfortunately, this means one can- 
not use the NIII or NV absorption profiles for constraining 
the ionization state or abundances of this system. At the 
same time, however, it is important to stress that these is- 
sues have little effect on the ionization fractions of C and 
Si as they are largely insensitive to the flux of photons 
with energies > 4 Ryd. 



5. SPECULATIONS AND FUTURE OBSERVATIONS 

The similarities and differences observed for the two 
metal-line systems toward J2233— 606 are compelling. The 
systems differ by only « 0.6 dex in iV(HI) and have com- 
parable total Hydrogen column densities. The inferred nu 
and ionization fractions are also in good agreement. Fur- 
thermore, the kinematic characteristics of the metal-line 
profiles closely resemble one another. The profiles are rela- 
tively narrow (< 100 km s _1 in velocity space) and exhibit 
only a modest asymmetry. One might expect, therefore, 
that while the z = 1.94 system may have a slightly lower 
AT(H) value, the two systems have a similar physical ori- 
gin. There is, however, a significant division between the 
two systems; the z = 1.94 system has a metallicity over 
20 times larger than the z — 1.92 system. We believe this 
striking comparison belies a fundamental difference in the 
physical nature of the two systems. 

To date, several models for the physical nature of the 
Lyman limit systems have been proposed. Several groups 
have demonstrated that a significant fraction of z f» 1 LL 

2 We adopt a lower limit to this ratio for the z = 1.92 system 
because we expect the NV A1238 profile may be an unidentified, 
coincident absorption line. 



Lyman limit systems to be associated with galactic or pro- 
togalactic systems. O n the other hand, numerical hydro- 
dynamic simulations ( Katz et al. 1996 ) indicate that LL 
systems may also arise in the large scale structure (LSS; 
e.g. filaments, walls) of the Universe. The z = 1.92 sys- 
tem, with its low metallicity, is a candidate for both of 
these models. If the system is associated with a galaxy, 
then either very little star formation (e.g. a dwarf galaxy) 
has taken place or the sightline intersects gas in the unen- 
riched outer parts of the galaxy. Also, the observation of 
two significant HI components would indicate substructure 
within a single galaxy, or perhaps that the sightline actu- 
ally intersects two galaxies. The fact that the metal-line 
profiles do not exactly coincide with either HI component, 
however, may be difficult to reconcile in terms of a galac- 
tic system where one might expect the HI gas to exactly 
trace the CII gas. At the same time, the observations 
can be explained within the framework of a LSS origin. 
This scenario, similar to the description of the Lya forest 
within N-body simulations, links the Lyman limit system 
with an unvirialized LSS in the universe, e.g., a filament. 
The metallicity is sufficiently small that the metals could 
have been injected at much earlier time by a galactic sys- 
tem as is proposed for the enrichment of the Lyq fores t 
( |Tytlcr et al. 1995| , |Cowie fc Songaila 199S| , |Lu et al. 1998) ) . 
The significant difference in the metallicity for the two HI 
components, however, argues for a local production of the 
metals (e.g. a hydrodynamic shock in the LSS). While it is 
not difficult to envision a LSS scenario which explains the 
two HI components and the misalignment of the metal- line 
profiles, more research on this topic is necessary to reveal 
the likelihood of these observations. Perhaps one will find 
it is a general prediction from the numerical simulations 
that one would not expect the HI gas to directly align with 
the metal-line profiles. 

In contrast, the high metallicity of the z — 1.94 sys- 
tem argues strongly against a LSS origin as it almost cer- 
tainly requires the presence of a significant star forming 
region. Because the system is within 30,000 km s _1 of 
the QSO, one must consider the possibility that it is an 
associated system. The fact that Ly-5 and the CIV pro- 
files are fully saturated indicates no partial covering and 
argues strongly against this interpretation. Furthermore, 
no NV absorption is evident as would be expected for an 
associated system. Therefore, we contend the absorber at 
z — 1.94 arises within a galactic system. In this scenario, 
the low N(H) and AT(HI) values indicate that the sightline 
has a large impact parameter. In fact, unless the system is 
the result of a very compact star forming region, we predict 
one will identify it with the STIS imaging and follow-up 
observations. If the responsible galaxy is identified, than 
our analysis affords a direct measurement of the physical 
properties of a z « 2 galaxy. 

By comparing our analysis with an examination of the 
imaging observations, one is able to investigate the the 
ISM in z w 2 galaxies (or in the case that no galaxy is 
identified, perhaps the properties of LSS). Unfortunately, 
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one cannot rely on photometric redshifts to distinguish be- 
tween z = 1.92 and z = 1.94, thus follow-up spectroscopy 
of all candidate galaxies will be essential. At the time 
of writing, we are aware of only one spectroscopic study 
flTrcssc et al. 1998[ ) of the STIS field. While their obser- 
vations identified no galaxies at z m 1.9, the observations 
were limited to I< 21. 
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REFERENCES 

Bergeron, J. & Boisse, P. 1991, A&A, 243, 344 

Boyle, B.J. 1997, AAO Newsletter, 83, 4 

Buries, S. and Tytler, D. 1998, ApJ, 507, 732 

Cowie, L. L. & Songaila, A. 1998, Nature, 394, 44 

Davidsen, A. F., Kriss, G. A., & Zheng, W. 1996, Nature, 380, 47 

Ferguson, H.C. et al. 1999, AJ, submitted 

Ferland, G. J. 1995, Hazy, a Brief Introduction to Cloudy, Univ. 

Kentucky, Phys. Dept. Int. Rep 
Gardner, J. P. et al. 1999, AJ, submitted 
Haardt, F. & Madau, P. 1996, ApJ, 461, 20 
Hamman, F. & Ferland, G. 1992, ApJ, 391, L53 

Katz, N., Weinberg, D.H., Hernquist, L., Miralda-Escudc, J. 1996, 
ApJ, 456, L57 

Lu, L., Sargent, W.L.W., Barlow, T.A., Churchill, C.W., & Vogt, S. 

1996, ApJS, 107, 475 
Lu, L., Sargent, W. L. W., Barlow, T. A. & Rauch, M. 1998, 

submitted to AJ 

Outram, P.J., Boyle, B.J., Carswell, R.F., Hewett, P.C., Williams, 

R.E., & Norris, R.P. 1998, MNRAS, submitted 
Pettini, M., Smith, L.J., King, D.L., & Hunstead, R.W. 1997, ApJ, 
486, 665 



Prochaska, J.X. 1999, ApJ, in press, istro-ph/9811357 
Prochaska, J. X. & Wolfe, A. M. 1997, ApJ, 486, 73 
Prochaska, J.X. & Wolfe, A.M. 1998, ApJ, 507, 113 
Prochaska, — J,X. & Wolfe, A.M., 1999, ApJ, in press, |astro 



ph/ 9,810301 

Savage, B. D. and Sembach, K. R. 1991, ApJ, 379, 245 
Savaglio, S. 1998, preprint 

Scaloy, K.M., Drinkwater, M.J., & Webb, J.K. 1998, ApJ, submitted 



astro-ph/980401S 
Steidel, CXT. 19 U U, ApJS, 74, 37 
Steidel, C.C. & Dickinson, M. 1992, ApJ, 394, 81 
Tresse, L., Dennefeld, M., .Pet it jean. P., Crjstiani, S., & White, S. 



1998, A&A, submitted, ( astro-ph/9812246;) 
Tytler, D., Fan, X.-M., Buries, S., Cottrcll, L, Davis, C, Kirkman, 

D., & Zuo, L. 1995, in QSO Absorption Lines, ed. G. Meylan 

(Springer- Ver lag), 289 
Valls-Gabaud, D. & Vernet J., MNRAS, in preparation 
Williams, R. et al. 1999, AJ, submitted 



